Hund's Rules for the N = Landau Levels of Trilayer Graphene 
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The N = Landau levels of ABC and ABA trilayer graphene both have approximate 12 fold 
degeneracies that are lifted by interactions to produce strong quantum Hall effects (QHE) at all 
integer filling factors between u = —6 and v = 6. We discuss similarities and differences between 
the strong-magnetic-field weak-disorder physics of the two trilayer cases, and between trilayer and 
bilayer cases. These differences can be used to identify the stacking order of high-quality trilayer 
samples by studying their quantum Hall effects. 
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I. INTRODUCTION 

Recent experimental advances^ have realized a 
new family of quasi two-dimensional electron systems 
(2DESs) consisting of two or more graphene layers. Many 
of the unusual properties of these systems are related to 
momentum-space textures in the low-energy quasiparti- 
cle bands, and to the correspondingly enhanced N = 
(zero energy) Landau level (LL) degeneracies that ap- 
pear when an external magnetic field is applied. As in 
the case of semiconductor 2DESs, Coulomb interactions 
break LL degeneracies when disorder is weak 2 -^ and fill- 
ing factors are close to integer values. This phenomenon 
is sometimes referred to as quantum Hall ferromagnetism 
because its simplest^ manifestation is spontaneous spin- 
polarization in the conduction bands of zincblende semi- 
conductors at odd integer filling factors. In few-layer 
graphene 2DESs, LL degeneracies are larger because of 
the presence of two low-energy valleys (K and K')§-—, 
and for N = also because of the momentum-space pseu- 
dospin textures. 

In bilayer graphene, the additional degrees of freedom 
can lead to canted antiferromagnetism^— , and valley 
and orbital order— ~— among other— ~— phenomena. The 
trilayer graphen o 19-22 system, which has become exper- 
imentally realized 2 ^— only recently, presents an even 
richer many-electron physics landscape. Both ABC (chi- 
ral) and ABA stacked trilayers are expected to have in- 
teger Hall plateaus at v = ±4(n + 3/2). The v = — 6 
and v = 6 plateaus are separated by a set of 12 approx- 
imately degenerate N = LLs, the trilayer duodectet, 
due to a combination of the 4-fold spin-valley degener- 
acy common to all levels, and trilayer momentum space 
pscudospin texture a 19 ' 20 . However, as we shall explain, 
the role of electron-electron interactions is quite different 
in the two cases because of differences ^ 26 ' 33 in the spatial 
and orbital contents of the N = LL wavefunctions. 

ABC trilayers have N — layer wavefunctions that 
arc localized on the A\ site in Fig. U for one valley 
and on the B3 site for the other. ABC trilayers arc 
like AB bilayers in that they have inversion symmetry 
protected 3 ^ degeneracy, but they differ in that they have 
N = LLs with n = 2 orbital character, in addition 
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FIG. 1: (Color online) Schematic electronic structure of few- 
layer graphene. The sites coupled by the strongest interlayer 
bonds characterized by the 71 hopping parameter are indi- 
cated by shading and have little weight in N — Landau 
levels. 



to n = and n = 1. Their quantized LL energies are 
~ B 3 / 2 ^n(n — l)(n — 2). This property causes LL 
mixing effects to vanish strongly at very high magnetic 
fields. The N = LLs of ABA trilayers, on the other 
hand, consist of n = 0, 1 orbitals on the layer symmetric 
states |Ai) + |A3) and I-B2), and only n = orbitals on the 
layer antisymmetric states \Ai) — \A 3 ) and \B\) — \B 3 ). 
The N LL spectrum can be viewed as a superpo- 
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- B 1 / 2 ^ and 

bilayer-like LLs with E n z> ~ B^/ n(n — 1). This structure 
leads to LL crossin g 22 ' 25 ' 35 and for N ^ increases the 
quantitative importance of LL mixing at stronger fields 
compared to the ABC case. 

Our description of the N = physics of ABC and ABA 
trilayer graphene starts from a model that includes only 
nearest neighbor intra-layer and inter-layer hopping and 
focuses on energies smaller than the inter-layer hopping 
strength 71. By treating other hoppings pcrturbatively 
at low-energies, it is possible to derive effective Hamilto- 
nians for N = LLs that are valid at realistic magnetic 
field strengths. In Section II of this article we discuss 
N = quantum Hall ferromagnetism in ABC trilayer 
graphene, predicting a simple set of Hund's rules that 
apply in the absence of inversion-symmetry breaking ex- 
ternal electric fields. We then address the phase tran- 
sitions that can be driven experimentally by applying 
electric fields across the trilayers using top and bottom 
gates. The interaction physics of the N = level in ABA 
graphene, which is discussed in Section III, is more com- 
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plex even in the absence of an external electric field. We 
conclude in Section IV with a discussion of prospects for 
interesting experimental discoveries in these systems. 



II. ABC-STACKED TRILAYER GRAPHENE 

The low-energy model for ABC trilayer o 19 ' 20 in the A\- 
B3 representation is: 
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In Eq.(TT|) 7r = hk + eA/c is the 2D kinetic momentum, 
7r = t z 7t x + in y where t z = ±1 represents valley K and 
K' , vq is the monolayer Dirac velocity, 7! ~ 0.5 e V 19 ' 27 , 
and 2ud is the site-energy difference between A\ and B3. 
This model is approximately valid for hvo /£b <C 71 where 



£b = 25.6 nm/y / B[T] is the magnetic length, i.e. for 
fields much smaller than ~ 400 T. For = 0, "Habc 
has three zero energy eigenstates for each value of t z . 
The corresponding eigenstates, which play the key role 
in ABC trilayer quantum Hall ferromagnetism have the 
form (0, </>„) for valley K and (<^> n ,0) for valley K' . Here 
n = 0, 1,2 is an orbital label and <j) n is the nth LL eigen- 
state of an ordinary 2DEG. The ABC trilayer duodectet 
is the direct product of this LL orbital triplet, a real spin 
doublet, and a which-layer (or equivalently which-valley) 
doublet. 

Gaps can open within the duodectet, giving rise to 
additional integer quantum Hall effects, due to single- 
particle terms that lift degeneracies or due to interac- 
tions that break symmetries. For example a Zeeman 
term splits spin levels by 2Ezm = <?«b-B = 0.116 x B[T] 
meV, while splits both layers and orbitals: E£ L = 
— r z Ud(l — n(fiw3/7i) 2 / 3 ). The notation used in this 
equation is motivated by viewing the ABC multilayer 
graphene Hamiltonian as the J = 3 generalization of 
Hamiltonians that apply to single-layer graphene for 
J = 1 and to bilayer graphene for J = 2. In each 
case hwj is proportional to the energy of the N = 1 
LL. hjji = (v / 2woft/7i^B) J 7i is the cyclotron frequency 
of quasiparticles in chirality J 2DES's. The Zeeman en- 
ergy is invariably small compared to the interaction scale 
e 2 /e^B, whereas Ud can be large enough to completely re- 
organize the duodectet. 

The full electron-electron interaction Hamiltonian in 
the duodectet subspace is readily available given the sim- 
ple form factors of the n = 0,1,2 Landau levels. In 
our discussion of N = quantum Hall ferromagnetism 
we will nevertheless use the Hartree-Fock approximation 
which is expected to be accurate 3 ^ at the integer filling 
factors of interest, and more physically transparent. We 
solve the Hartree-Fock equations in the duodectet sub- 
space allowing any symmetry, other than translational 
invariance, to be broken. The Hartree-Fock Hamiltonian 
is expressed in terms of the duodectet density-matrix 



which must be determined self-consistcntly at each filling 
factor. When we include only interactions 3 - 7 , within the 
N = manifold we find the quasiparticle Hamiltonian 
with Hartree and exchange interaction contributions is 

(cms\Hfi c \Pn's') = (E£ L 5 SS , - E zu <t z ss >) <W<W 

: A ,9niS 'r> , i2\ 



where we have used the labels n — 0,1,2 for orbitals, 
index s for spin, and the Greek indices a,/3 = B{K) 
or T(K') for layer(vallcy). Here Ep = e 2 /e(.B charac- 
terizes the overall strength of the exchange interactions, 
-Eh = (2g?/2£b)-Ef captures the capacitive electrostatic 
energy associated with shifting charge between the top 
and bottom layers, and d = 0.335 nm is the separation 
between adjacent graphene layers. The layer projected 
charge density is 



/ . "an 



and the density matrix is 
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The dimensionless exchange integrals X"£„.„, reflect 
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where ry a/3 becomes 1 for a = /3 and e _2kd for a / /). 
Here F nn i(k) are the ordinary^ 2DEG LL form factors 
which can be expressed in terms of associated Laguerre 
polynomials that capture the spatial profile of the LL 
wavefunctions. The exchange integrals are non-zero only 
for n — n' = 712 — ni; numerical values of the relevant 
form factors for the case d — arc listed in Table fl] 

The self-consistent solutions of the N = ABC 
Hartree-Fock equations for balanced (ud — 0) ABC tri- 
layers are summarized in Fig(2] for a typical field of 
B = 20 T. Gaps appear at all intermediate integer fill- 
ing factors, and as expected are always much larger than 
Ezm- The electronic ground state at duodectet filling- 
factors from v = —6 to 6 follows a Hund's rule behav- 
ior progression in which spin polarization is maximized 
first, then layer (valley) polarization to the greatest ex- 
tent possible, and finally LL orbital polarization to the 
extent allowed by the first two rules. In these calcula- 
tions, the quasiparticle energies of spins along the field 
direction (say spin f) are much lower than those for spin 
J. states because the spin-splitting is exchange enhanced. 
Since d/ £-q is always small in trilayer graphene systems, 
the exchange integrals that favor spin-polarization are 
only slightly larger than those that favor layer polariza- 
tion. Because these differences and Zeeman energies are 
both small, inter-valley exchange 3 - 9 - (which is neglected 
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TABLE I: Summary of the numerical values of exchange integrals. The intralayer integrals are independent of magnetic fields 
while the interlayer ones decline slowly as the field strength increases. Numerical values are reported for B — 20 T and B = 30 
T. We have defined indices a = |n — n-2 1 , b — min(m, n') and c = min(n, 722). Exchange integrals are symmetric between b and 
c. For the interlayer exchange integrals the values reported here are for neighboring layer separation d = 0.335 nm. 
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in the continuum model we employ here) can potcntialy 
alter the competition between spin and layer polarization 
compared to the results reported in FigJ5] We return to 
this point in the discussion section. For balanced ABC 
trilayers, layer-order appears as spontaneous inter-layer 
phase coherence^ - — rather than layer polarization, and 
to Hartrcc-Fock quasiparticles that have symmetric (S) 
or anti-symmetric (AS) layer states. The three LL or- 
bitals are filled in the order of increasing n because ex- 
change integrals decrease with n as indicated in Table |U 
It follows that the first six filled LLs are \S, t), \S, 1 1) 
and \S, 2 f) followed by the AS counterparts; the next six 
filled LLs are the spin i states in the same order. Here 
S and AS are orthogonal coherent bilayer states, chosen 
as symmetric (S) and antisymmetric (AS) states as a 
matter of convenience. 

These Hund's rules imply that the interaction driven 
integer quantum Hall states are spin and pseudospin po- 
larized at all the intermediate integer filling factors from 
v = —5 to 5, as summarized in FigJ5] All 11 states are 
spin f polarized fcrromagncts with a maximum spin po- 
larization at v = 0. The LL orbital dependence of the 
microscopic Hamiltonian leads to orbital polarization ex- 
cept at filling factors that are multiples of 3. In the ab- 
sence of an interlayer potential which breaks inversion 
symmetry, the layer pseudospin symmetry is broken by 
states with spontaneous interlayer coherence which open 
a gap at no cost in Hartree energy. The Hartree en- 
ergy and the difference between intralayer and interlayer 
exchange are both the same order ~ d/l-Q ■ e 2 /£-B, but 
it turns out that Hartree energy slightly dominates over 
exchange difference favoring spontaneous coherence state 
over competing layer polarized states. 

We have so far neglected the strong dependence of 
duodectet quantum Hall fcrromagnets on an interlayer 
electric field. A very small potential difference between 
the top and bottom layers is sufficient to destabilize 
spontaneous-coherence states in favor of layer polarized 
states. Quantum phase transitions vs. u,i are expected 
at all filling factors. For v = 0, for example, there is a 
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FIG. 2: (Color online) Upper panel: filling factor dependence 
of the Hartree-Fock energies of occupied and unoccupied LLs 
for a balanced ABC trilayer at 20 T. Energies are in units of 
(n/2) 1 ^ 2 e 2 /eIb- Unoccupied LLs are sextets (blue), triplets 
(green), or singlets (cyan); occupied LLs are doublets (ma- 
genta) or singlets (red); LLs at v = ±6 are split only by the 
small Zeeman splitting. These energies include only exchange 
interactions within the N = level. Particle-hole symmetry 
within the N — space is recovered if exchange interactions 
with remote Landau levels are included; these interactions 
do not influence the energy spacings between Landau levels. 
Lower panel: filling factor dependence of the polarizations of 
up spin (blue) and layer (red), and population differences be- 
tween n = 2 and n — (green) and n — 1 (cyan) LL orbitals. 



first-order phase transition from the spin polarized ferro- 
magnetic state to a layer(valley) polarized state when the 
potential difference between the top and bottom layers 
becomes dominant over the exchange induced spin split- 
ting. For B = 20 T (e — 1) we find that the layer polar- 
ized state has the lowest energy if the potential difference 
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between the top and bottom layers exceeds 0.2 eV. The 
gap between the highest occupied LL \T, 2,f) and the 
lowest unoccupied LL \B, 0,1) is 0.46 eV for Ud = and 
is reduced to 0.28 eV at the phase transition. We find 
that the critical screened electric field for this first order 
phase transition is ~ 15 mV/(nm-T) for e = 1, compara- 
ble to the experimental value of the corresponding tran- 
sition field in bilayers42r— The linear dependence of the 
critical field on B follows from the fact that all relevant 
microscopic energy sales (the bias supported energies En 
and Ell j the difference between intralayer and interlayer 
exchanges, and the Zeeman energy £zm) are all linear 
functions of B. 

In high mobility and low disorder suspended sam- 
ples, we anticipate that the v = ±6, states and the 
bias supported v = ±3 states, persist down to nearly 
zero magnetic field exhibiting spontaneous quantum Hall 
effectsj&^i ABC trilayers are highly susceptible to broken 
symmetries that are accompanied by large momentum 
space Berry curvatures and different types of topological 
order4 Remote weaker hopping processes^ in ABC tri- 
layers that we have ignored so far could influence the field 
to which the quantum Hall states survive. The most im- 
portant of these appears likely to be the trigonal warping 
energy scale^ 7 meVJ^ The study of weak field quantum 
Hall effects in ABC trilayers should shed light on the na- 
ture of the B = ground stated 



III. ABA-STACKED TRILAYER GRAPHENE 

Unlike ABC trilayer quasiparticles with pure chirality, 
the band structure of an ABA trilayer consists of a mass- 
less monolayer (J = 1) and a massive bilayer (J = 2) 
subbands. Unbiased ABA trilayers have mirror symme- 
try with respect to the middle layer and their low energy 
physics is governed by2£ 
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J=2 



where the J = 1 subbands are layer antisymmetric states 
\Ax) - \A 3 ) and \Bi) - \B 3 ) while the J = 2 subbands 
are layer symmetric states |Ai) + \A 3 ) and \B%). Note 
that the effective interlayer hopping energy is enhanced 
by a factor of \/2 compared to the AB bilayer case. 
Again, the mirror symmetry leads to spontaneous coher- 
ence states. The next-nearest-neighbor tunneling pro- 
cesses play a more essential role than in the ABC tri- 
layer case. Here 72 = —20 meV, 75 = 40 meV and 
5' = 50 me V 22 ' 25 lead to band gaps for the J = 1 and 
J = 2 chiral branches separately, but no direct gap over- 
all. When Zeeman splittings and interactions ignored, 
the duodectet has four-fold degenerate J = 2 eigenener- 
gies at e ~ — 10,0meV and two- fold degenerate J = 1 
eigenenergies at e ~ 10, 30meV. These energies are equal 
to the B = band energy extrema and are independent 



of magnetic field. The single-particle splitting of N = 
can therefore lead to quantum Hall effects at v = —2, 
v = 2 and v = 4. 

The self-consistent Hartree-Fock Hamiltonian that de- 
scribes the broken symmetry ABA trilayer N = 
duodectet is 

(i Hi s\V.%A\j n 3 s') = (E LL <$ SS ' - ^ZMCrfsO ^n.n^ij 

, Eh 



-A B2 (2<5b 2 ,i - l)<W<Wi J £ij 

^E^n 2 ,n 1 n J Ai^'^, (7) 



where the LL index ri; range depends on the atomic 
orbital i, i.e., n% = only for J = 1 branch while 
rii = 0, 1 for J = 2 branch. The exchange integrals 
X 1 ' 3 are still defined by Eq.© but with a more gen- 
eral definition if 3 = X]m c m^m/K)j where Vb,i,2 respec- 
tively denote intralayer, nearest interlayer, and next- 
nearest interlayer Coulomb interactions, and c m is ob- 
tained from the interaction matrix element decomposi- 
tion (ij\V\ij) = X)m c m^m- Otherwise the notation in 
Eq.© is the same as in the ABC case. Since it is already 
quite complicated, we focus on the balanced cas o 47 i 48 
with Ud = 0. -Ell is therefore absent and replaced by 
E LL which is given by the tt — eigenvalues of Eq.©. 

We find that the ABA duodectet intermediate filling 
factor states follow the Hunds rule for maximum spin- 
polarization, but are otherwise strongly influenced by 
the single-particle gaps. The J = 2 LLs, which have 
lower single-particle energies are filled first. When the 
magnetic field is smaller than B c \ = 17 T (s = 1 here- 
after), the first six filled LLs are all majority spin levels: 
|Ai + A 3 t) with n = and 1, then \B 2 t) with n = 
and 1, and lastly \A\ — A3 f) and \B\ - B 3 f) with 
n = 0. The next six filled LLs are the spin i states 
in the same order. Within a given spin the LL filling 
pattern follows the order suggested by the single-particle 
physics. The Hartree-Fock theory predictions for field 
strength B — 10 Tesla are summarized in Fig(3]Ja). At 
v = —5 and v = —4 f) LLs are occupied first. In 

addition to being favored by the single-particle energies, 
this choice gains exchange energy and minimizes Hartrcc 
energy. At v — —3 the I-B2 0, i") LL is then occupied 
in agreement with the single-particle Hamiltonian, but 
also lowering the interaction energy. Notice the the gaps 
between occupied and empty states, which are impor- 
tant for transport, arc greatly enhanced relative to the 
single-particle gaps at all intermediate filling factors. The 
smallest Fermi level gaps occur at filling factors v = — 5 
and v = —3, for which occupied and empty states are 
separated by a gap between n = and n = 1 orbitals 
with the same layer structure. These two filling factors 
should therefore have the weakest quantum Hall effects 
in the weak-disorder samples to which the present calcu- 
lations apply. 

Although the LL filling order is universal for v > — 3, 
there is some field dependence for v = —5, —4 and —3 
as indicated in FigJ3] When the field strength exceeds 
B c i = 17 T a LL reordering occurs at v = —4. We find 
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FIG. 3: (Color online) (a) Filling factor dependence of the 
Hartree-Fock (HF) energies (eV) of occupied and unoccupied 
LLs for a balanced ABA trilayer at 10 T. Unoccupied LLs are 
quartets (green), doublets (blue), or singlets (cyan); occupied 
LLs are doublets (magenta) or singlets (red), (b)-(d) Field 
dependence of the HF energies (eV) of the lowest four LLs at 
v = —5, —4 and —3. Results for magnetic field B equal to 
15 T are shown in (b), for 25 T in (c) and for 35 T in (d). 
Magenta and red denote the n = and 1 orbitals respectively 
of LL \A\ + Ai f); green and blue denote the n = and 1 
orbitals respectively of LL |i?2 t)- The Zeeman splitting is 
small enough to ignore in this figure and we assumed that 

£ = 1. 



that the pair of LLs \B 2 t) are then filled before the 
\A\ + A^ t) pair. This transition occurs when Coulomb 
interaction physics overwhelms the single-particle next- 
nearest layer tunneling effect. More localized layer or- 
bitals are preferred when interactions dominate because 
intralayer exchange energy gains dominate interlayer ex- 
change and Hartrec energy costs. When the field is fur- 
ther strengthened to B c2 = 33 T LL \B 2 t) is filled 
before | A\ + A3 t) ; at v = —5. The two critical fields 
are roughly determined by the condition -Eh/2 ~ AE[^ h 
which coincides approximately with the self-consistent 
numerical results for B c \ and B c2 . 

We have focused here on interaction physics in the 
N = LL's of ABC and ABA trilayer graphenc. Interest- 
ing effects are nevertheless particularly likely for larger N 
in the ABA case. Unlike ABC trilayers where only J = 3 
bands are present at low energies, both J = 1 and 2 sub- 
bands appear at all energies in ABA trilayers, leading to 
LL crossin g 2 2 i 25 ~ between t he two chiral branches when 
u)i/u>2 = \Jn 2 {n 2 — l)/2m is satisfied. We anticipate 
that Coulomb physics will open gaps at the LL crossing 
points, whose characteristics would depend on the orbital 
indices of the crossing LLs. Because of the LL crossing 
effects present even without interactions, the sequence of 
the plateaus in ABA trilayers highly depend on the mag- 
netic field strength. 



IV. DISCUSSION 

Single-layer, bilayer, and ABC and ABA trilayer 
graphene are a family of related but distinct two- 
dimensional electron systems. Single-layer and bilayer 
graphene have four and eight fold Landau level degenera- 
cies that are lifted by interactions in the quantum Hall 
regime when disorder is sufficiently weak, giving rise to 
a rich variety of broken symmetry states. In this pa- 
per we have explored the way in which interactions alter 
the quantum Hall effect in both ABC and ABA trilayer 
graphene in the range of filling factors between v = — 6 
and v = 6, over which their twelve-fold N = Landau 
level is filled. As in the single-layer and bilayer cases our 
calculations suggest that quantum Hall effects will even- 
tually occur at all intermediate filling factors as sample 
quality improves to lessen the role of disorder. There are 
a number of clear differences between the ABC and ABA 
cases which could provide a practical alternative to Ra- 
man scattering for identifying the stacking orders of high 
quality trilayer samples: i) Because the N — layer or- 
bitals are localized either in the top or bottom layer in the 
ABC case, while they are often spread across two or more 
layers in the ABA case, interaction effects are stronger for 
ABC trilayers and more likely to induce intermediate fill- 
ing factor quantum Hall effects, ii) Perpendicular electric 
fields have relatively simple consequences in ABC trilay- 
ers, opening up gaps between layers that will strengthen 
the quantum Hall effect at v = ±3 and v = 0, but have a 
much more complex behavior with many quantum phase 
transitions in the ABC case. For ABC trilayers there is 
a single sharp phase transition from a spin-polarized (at 
high field) or antiferromagnetic (at low field) v = state 
to a layer polarized state with a critical field strength 
that is linear in B. iii) In ABA samples on substrates, 
the single-particle splitting of N = LLs might lead 
to quantum Hall effects at 1/ = - 2, v = 2 and v = 4 
with gaps ~ 10 — 20 meV before the interaction driven 
quantum Hall states appear at larger fields, iv) There 
are phase transitions in ABA trilayers at v = —5, —4, —3 
between low-field states in which single-particle physics 
determines the order in which LL's are occupied, to high- 
field states in which more localized orbitals with larger 
exchange energies are occupied first, v) In high quality 
suspended ABC trilayers, v = 0, ±3, ±6 states persist 
down to zero magnetic field and smoothly connect with 
the electron-electron interaction driven broken symme- 
try states^. We expect that these differences will become 
apparent in future experiments. 

In this paper we have reported on Hartrec- Fock ground 
states and quasiparticle energies as a function of filling 
factor predicted by calculations which were restricted 
only by the assumption of translational invariance but 
otherwise allowed any symmetries of the Hamiltonian to 
be broken. We find that in the strong-magnetic-field 
weak-disorder case the LL filling order tends to follow 
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a pattern in which states of a particular spin arc filled 
first, and within a spin states with different valley content 
next. For ABA trilaycrs, however, single-particle physics 
which places a bilayer-like group of states lower in en- 
ergy than a monolayer-like group of states, can play an 
important role. The rule which favors spin-polarization 
over valley (or equivalently layer) polarization is weak be- 
cause the layers in this system are quite closely spaced, 
and it is possible that this rule will be overturned by lat- 
tice effects - as it appears to be in bilayers^. When M 
of the 12 N — Landau levels are occupied, there will 
be a set of M x (12 — M) collective modes whose ener- 
gies approach differences between quasiparticle energies 
in the limit of long-wavelength. Many of these collec- 
tive modes, which can be studied^ using time-dependent 
versions of the calculations reported in this paper, will 



couple to long-wavelength electromagnetic probes which 
would provide information about the many-body ground 
state that is complementary to the transport activation 
gaps most related to the present study. 
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